Attosecond light pulse induced photo-association 
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We explore stimulated photo-association in the context of attosecond pump-probe schemes of 
atomic matter. An attosecond pulse - the probe - is used to induce photo-association of an electronic 
wave packet which had been created before, typically with an attosecond pump pulse at an atomic 
center different from the one of photo-association. We will show that the electron absorption is 
maximal for a certain delay between the pulses. Two ways of enhancing and controlling stimulated 
photo-association are proposed, namely using an additional infrared pulse to steer the electronic 
wave packet and using a train of attosecond pulses instead of a single pair. A direct application of 
ultrafast stimulated photo-association is the measurement of atomic distances. 

PACS numbers: 32.80.Rm, 42.50.Hz 



I. INTRODUCTION 



Attosecond laser pulses hold the promise to be a tool 
to image structure of matter on the atomic scale with 
unprecedented temporal and spatial resolution [H, 0, 9 • 
Nowadays, attosecond pulses with a typical duration of 
electronic motion in atoms or molecules [H can be pro- 
duced and characterized in a controlled way Q , and are 
therefore ready to be used for mapping and controlling 
electronic motion [|| 0] ■ First applications of attosecond 
pulses have focused on the generation of attosecond elec- 
tron wavepackets, e.g., for real time observation of tun- 
neling [H], the interference of wavepackets using attosec- 
ond trains @, or to obtain autoionization widths [icj |. 
They have also been used for imaging molecular orbitals 
[Til ] . All the applications show that the role of attosec- 
ond pulses has been restricted so far to study ionization 
or interference processes. Here we propose a new applica- 
tion of attosecond pulses: the recapture of a continuum 
electron by an ion. induced by an attosecond pulse. This 
constitutes a new process of photo-association, namely 
induced photo-association (IPA) . Although it is a process 
of stimulated photo-absorption triggered by the attosec- 
ond photon field, it has low probability due to the small 
timespan of the attosecond pulse during which the IPA 
is possible. 

Hence, after the demonstration of the basic process of 
IPA we will also discuss how the reabsorption probabil- 
ity can be enhanced. This can be achieved by using an 
attosecond pulse train, which delivers a series of pairs of 
pump-probe pulses. The drawback of such a technique is 
that the time delay between two attosecond pulses is typ- 
ically not variable in an experiment. However, adding an 
infrared (IR) pulse can achieve in the context of IPA the 
same effect as changing the delay between the attosecond 
pulses, namely controlling the time when the electronic 
wavepacket arrives at the atomic center by which it is 
supposed to be captured. This will be shown explicitely 
for one attosecond pulse per IR cycle, which is experi- 
mentally achievable 

Since this work is mainly intended to introduce the 



concept of IPA, we restrict ourselves to the simplest pos- 
sible scenario, two ions and an electron, and to be spe- 
cific, we take ■ Furthermore, we carry out fully time- 
dependent quantum calculations in one and two spatial 
dimensions. 

The paper is organized as follows. In section II we in- 
troduce our Hamiltonian and describe how we solve the 
electron dynamics. In section III we describe the basic 
IPA process with two attosecond pulses. In section IV we 
discuss how the addition of an IR pulse changes the dy- 
namics. Section V presents the results for an attosecond 
pulse train in combination with an IR pulse, and sec- 
tion VI contains conclusions and outlook of the paper. 
Atomic units are used unless stated otherwise. 



II. TIME-DEPENDENT ONE-ELECTRON 
QUANTUM DYNAMICS IN A LANDSCAPE OF 
ATTRACTIVE POTENTIALS 

As mentioned before, for the sake of clarity of the con- 
cept, we consider the simplest possible dynamical situa- 
tion for IPA to occur in a pump-probe setting with the 
ions (protons) fixed in space. We have in mind the ide- 
alized situation that an electron, bound to one proton 
(forming hydrogen) is ionized by the first atto pulse. The 
electron packet propagates and part of it is recaptured 
subsequently at another proton, assisted by a second atto 
pulse. The recapture is maximal if the probe pulse comes 
at the time when the center of the wavepacket is at the 
second proton. 

The fixed nuclei assumption is a twofold but well jus- 
tified approximation. Firstly, the situation of a hydro- 
gen atom and a proton separated by 40 a.u. as con- 
sidered below, is naturally created by exciting the H£ 
molecule from its ground state, with equilibrium inter- 
nuclear distance i?o = 2, to a dissociative state. This 
leads to a maximum nuclear kinetic energy determined 
by the energy of this state, which is V2 pCT „ (Ro) = —0.169. 
At a distance R = 40, the potential energy approaches 
V2 p(Tu (oo) = —0.5. The difference determines the max- 
imum kinetic energy, Mv 2 = V2 p(Tu (i?o) — V2 p(Tu (oo). 
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Hence the relative velocity of the two nuclei with mass 
M = 1836 is v ~ 0.0134, so that the nuclei move less than 
half an atomic unit during a typical pump-probe inter- 
val of about 50 a.u.. Even smaller is the motion due to 
the IPA process itself, i.e., due to the Coulomb repulsion 
of the two protons after the first (ionizing) attosecond 
pulse. 

We assume the light field to be linearly polarized along 
the x-axis (unit vector x), 



E(t) = xE 
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(1) 

where t' = t — At and At is the delay between the two 
attosecond pulses, which have a central frequency ui and 
a duration of N cycles. The Hamiltonian reads 

P 2 

H=fj + V+p x A(t)/c=H + Px A(t)/c, (2) 

where A(t) = — c J* E(t') dt' . Furthermore, we consider 
only the case where the two protons are fixed on the x- 
axis at distances ±Xq with xq = 20, so that the soft-core 
potential reads in two dimension V(x, y) = W-(x, y) + 
W+(x,y) with 



W±(x, y) = -[(xT x ) 2 +y 2 + e] 



(3) 



The softening parameter e, adjusted to provide a ground 
state energy of E g = —0.5 for the isolated H atom, is 
e\D = 2.0 and €2D — 0.63 in one and two dimensions, 
respectively. The ground state energy corresponds to an 
isolated hydrogen atom, which is approximately true for 
the second proton at a distance of 40 atomic units. To 
localize the state at one nucleus, one has to take super- 
positions of the ground states of gerade and ungerade 
symmetry, ipg/u i which are formally the ground and first 
excited state of the system. This is well known from 
ion- atom collisions, 



(4) 



where we define ip+ to be localized at the right well 
W+(x,y) at x = +20 of Eq. © 

The initial state is propagated numerically in time un- 
der the Hamiltonian of Eq. @ using the Crank-Nicholson 
method with a time step of St = 0.05. The grid ranges 
from —100 to +100, with a spatial step Sx = 0.1. The 
cigenstatcs ip g / u from Eq. ^ are computed using imag- 
inary time-propagation under the Hamiltonian Hq from 
Eq. ©. 



III. TWO ATTOSECOND PULSES 

In this section we will investigate under which condi- 
tions induced photo-association can be realized within 
the framework of a dual attosecond pulse pump-probe 
scenario [131 ]. For simplicity, we will assume that at the 
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FIG. 1: Electron density as a function of time, under the ac- 
tion of a single attosecond pulse of 40, 60 and 80 nm central 
wavelength and N = 6 cycles (the central time of the pulse is 
approximately 16.5, 25 and 33 a.u., respectively). The wells 
with potentials W± are located at ±20 a.u., marked by hori- 
zontal dashed lines. The electronic ionization path predicted 
according to Eq. (((J (solid) and Eq. (0 (dashed) is indicated. 



beginning the electron is always well localized around 
xq in the well W+ with initial wavefunction ip+(t = 0). 
Asymmetric dissociative states can be realized (with an 
electron localization probability up to 84%) by means 
of laser-assisted dissociation Q. We use three different 
wavelengths (40 nm, 60 nm and 80 nm), and high in- 
tensity attosecond pulses (7 = 10 16 W/cm 2 ) in order to 
obtain good statistics for IPA. 

For optimum recombination, the main velocity v of the 
electron wavepacket, the delay At between the pulses and 
the distance 2xq between the two atomic centers must 
fulfill roughly the classical relation v = 2xq/ At. More 
precisely, since the electron absorbs one photon of energy 
lu from the pump pulse, energy conservation can be used 
to determine the electron velocity of the center of the 
wavepacket classically for a calculated electron trajectory 
r(t). With the Hamiltonian H from Eq. ^ we have, for 
t > 2-kN/lu, 



P 2 (t) 



+ V(r(t)) =<J + Ei 



Hence, 



v(t)=p(t) = {2[ L o + E l -V(r(t))]} 1/ \ 



(5) 



(6) 



and ignoring the potential energy for electrons with large 
energy in the continuum, the velocity becomes time- 
independent, 



vit) 



,(«) = 



= [2(cJ + £y] 



1/2 



(7) 



Results for the action of one attosecond pulse with A =40 
nm (Hlj =1.14), 60 nm (0.76) and 80 nm (0.57) are shown 
in Fig. [T] for Ei = E g . The attosecond pulse has a dura- 
tion of N — 6 cycles, and an intensity I = 10 16 W/cm 2 . 
The classically calculated path of the electron is shown 
for both, free electronic motion (dashed lines, Eq. ([7])) 
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and taking into account the nuclear potential (full lines, 
Eq. ©). While the free electron approximation is rea- 
sonably good for the most energetic electrons (40 nm) 
with excess energy of u> + E g = 0.614 (at a wavelength 
of 40 nm), it is not valid any more for the lower excess 
energy of 0.045 (80 nm). 

The time evolution of the eigenstates under the influ- 
ence of the total Hamiltonian given by Eq. ^ is 



0.025 



1>±(r,t) = U(t)1>±{r,0), 



(8) 



where U{t) is the time-evolution operator. We can also 
define the time evolution of the eigenstates under Hq, 
which is time-independent, 



iMr,*) = U o m±(r,0) = cxp-" ' V±(r,0). (9) 

As a measure for the fraction of the electron wavepacket 
which gets recaptured at the left well we calculate 
the overlap of its eigenfunction ip_ (x, t) (see Eq. ((4])) with 
the time-dependent wavepacket 



c_(f) = {i>_(r,t)\ii>(T,t)), 



(10) 



which is formally a cross-correlation function. 

The probability for induced photo-association (IPA) in 
the left well is therefore given by 



IPA 



(11) 



The final population of W~ is P[pA(tf), where tf is the 
time at the end of the second attosecond pulse, which is 
delayed by At with respect to the first one. As can be 
seen from Fig. [2 an optimum delay At* for maximum 
IPA exists. At* is reasonably predicted by the high pho- 
ton energy approximation for 40 nm, but not for lower 
energies (dotted lines). However, the inclusion of the nu- 
clear potential (Eq. ([6]), dashed lines) provides a good 
prediction for all three energies, with differences with re- 
spect to the real value of around 8% for 40 nm, and 10% 
for 60 and 80 nm. Note, that larger wavelengths pro- 
vide a worse accuracy in the predicted optimum delays, 
but correspond to higher absorption probabilities, which 
might be useful experimentally. 

The experimental observation of IPA is possible by de- 
tecting the outgoing proton or electron flux. For optimal 
time delay sufficient absorption at W- will lead to a de- 
crease in the number of protons as well as the number of 
electrons detected in this direction (negative x-axis). 

To study the asymmetry in the electronic emission, we 
use the continuity equation 



a|*(r,*)F 

dt 



-Vi(r,t), 



(12) 



which relates the change in the electronic density as a 
function of time with the particle flux 



j(r,t) = »( **(r,t) — V*(r,t) 
im 



(13) 




50 100 
Time delay (a.u.) 

FIG. 2: Probability of the induced photo-association at W- 
as a function of the delay between the pulses, for three photon 
wavelengths: A =40, 60 and 80 nm (one-dimensional calcu- 
lations). The predicted optimal delay according to Eq. © 
and Eq. Q is indicated with dashed and dotted lines, respec- 
tively. The arrows highlight the difference between the real 
and calculated values if the potential energy is ignored. 



In one dimension, integration of Eq. (JT2J) in space leads 
to 



d\^(x,t)\ 2 , dN . T , , T . 
1 y m n dx = — =j(-L,t)-j(L,t), (14) 



where N is the number of electrons. Using Eq. 
for the right hand side of Eq. (fT4"|) yields 



dN 
~dt 



ON, <9AL 



at 



dt 



again 



(15) 



providing the ionization rate ^ in terms of particles 
N + escaping to the right passing the position +L and 
the ones N- escaping to the left passing —L. We define 
the parameter 



A 



AN + - AN- 
AN+ + AN- 



(16) 



for asymmetric escape. If the electron is originally lo- 
cated in the well W+ (at +ro), the flux should ideally 
be symmetric if there is no absorption at VF_ . With ab- 
sorption at W- located at — ro, there should be less flux 
through — L than through L, so that A > 0. In Fig.[3]the 
asymmetry parameter for the two wells system is shown 
(full line) together with the IPA probability (dashed line, 
see also Fig. [3]) for comparison. Indeed, A has a maxi- 
mum if the population at the second well is maximal. 

Note, that the asymmetry parameter is never zero. 
The second well reflects part of the ionized wavefunction, 
resulting in a higher value of ionization flux towards the 
side where W+ is located (positive x-axis). For this rea- 
son the experimental detection of the reabsorption in W- 
should be based on the shape of the asymmetry param- 
eter, and not on its absolute value. 

Our one-dimensional calculations suggest that electron 
population can be efficiently transferred to the second 
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FIG. 3: Full line: Asymmetry parameter as a function of the 
time delay (Eq. (|16[) ). for A = 80 nm and the same laser 
parameters as in Fig. [2] Dashed line: IPA probability (popu- 
lation at W-), shown for comparison. 
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FIG. 4: The same as Fig. [2] but in two dimensions. 



well W- by optimizing the time delay between two at- 
tosecond pulses. The induced photo association probabil- 
ity as a function of the time delay can be used to measure 
the distance of the nuclei in the molecule. Furthermore, 
we have seen that the left and right asymmetry in the 
ionization signal could be a natural way to observe these 
effects in an experiment. 

To substantiate these findings we have carried out 
more realistic 2D calculations, where spreading of the 
electronic wavepacket in the coordinate transversal to 
the laser polarization is naturally included. Results 
shown in Fig. 2] are fully analogous to those of the one- 
dimensional case, although the IPA probabilities are ob- 
viously smaller, since the outgoing electron has now ac- 
cess to the transverse direction. The differences between 
obtained and predicted maxima according to Eq. ^ are 
now 12% for 80 nm. 6% for 60 nm and 9% for 40 nm. 



IV. TWO ATTOSECOND PULSES PLUS 
INFRARED FIELD 

So far we have demonstrated that the time delay be- 
tween two attosecond pulses can be used to measure 
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FIG. 5: Probability of the induced photo-association at W- as 
a function of the carrier-envelope phase of the infrared pulse. 
The insets show the total field for = and n. 



atomic distances, or to optimize the electronic transfer 
between two nuclei. However, with the current state-of- 
the-art technology, it is still difficult to isolate two pulses 
and choose the time delay between them. 

Attosecond pulses are created [l[ through High Har- 
monic Generation produced by the interaction of an IR 
laser field with an atom or molecule. With this technique 
one can easily generate Attosecond Pulse Trains (APT) , 
with one or two [ljj attosecond pulses per infrared cycle 
of the generating field. For a Ti:Sa laser (A =800 nm), the 
time delay between attosecond pulses is fixed to T ~110 
or T/2 ~55, and therefore the technique described in the 
previous section cannot be directly applied. 

To overcome this limitation we will introduce, in ad- 
dition to a fixed and realistic delay between attosecond 
pulses, a synchronized IR field. Due to the way the APT 
is generated, a synchronized IR pulse is always available. 

The phase of the IR laser pulse can be chosen to either 
accelerate or decelerate the electron when it is released 
at the first well W+ [HI, EH and therefore it will influ- 
ence the temporal and spatial overlap between the ion- 
ized electron and the second well W- at the time of the 
second attosecond pulse. 

In other words, an appropriate infrared field, applied 
while the electron travels the distance between the nuclei, 
can alter the electron trajectory such that the electron 
arrives at W- when the second attosecond pulse is ap- 
plied, thereby optimizing IPA. 

Our electric field has from now on three contributions: 
two Gaussian attosecond pulses of intensity Ijjv = 10 15 
W/cm 2 and FWHM = 2T, and an infrared field 



Em(t) = Em sin(wiR,7j + </>) sin' 



/ bj IR t \ 
\2N m ) 



(17) 



with N m = 8 and I m = 10 13 W/cm 2 . We have re- 
duced the UV intensity by a factor of 10 compared to 
the previous section, to work out more clearly the IPA 
enhancement produced by the IR field. The two attosec- 
ond pulses, separated by an infrared period At; = Tjr, 
are located symmetrically with respect to the center of 
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FIG. 6: Bending of the electron trajectory caused by the 
infrared field (A =800 nm, I = 10 13 W/cm 2 ). Left: with 
carrier-envelope phase <f> = in the infrared pulse (see insets 
in Fig. [5}. Right: the same situation when (f> = ■k. Dashed 
black line: position of the nucleus W-. Dashed white line: 
position of the center of the second attosecond pulse. Full 
line: expected trajectory following Eq. (|6]). 



the infrared. The central wavelength of the attosecond 
pulses is A = 65 nm (w =0.7). This implies that the op- 
timal delay (following Eq. |(7J)) would be t = 2xo/v ~ 68, 
much less than the actual delay of T = 110 between the 
pulses. Therefore, we expect enhanced recombination at 
W- if the IR field decelerates the outgoing electron. To 
this end we vary the carrier-envelope phase of the in- 
frared pulse, <fi in Eq. (fT7|) . This changes the value of 
the IR field at the moment when the electron is ionized 
by the first attosecond pulse and is equivalent to chang- 
ing the time at which the attosecond pulses occur during 
the infrared pulse. The latter has already been achieved 
experimentally (l6| . 

For the present intensity, the infrared pulse by itself 
causes negligible absorption at W- of only Pipa ~10 -13 . 
This ensures that IPA comes from the attosecond pulse 
only, while the infrared pulse has exclusively the role of 
accelerating or decelerating the ionized electrons. Results 
for IPA in the second well as a function of the carrier- 
envelope phase of the infrared pulse are shown in Fig. [5] 
The insets show the position of the attosecond pulses 
with respect to the infrared pulse, for the cases with 
4> = and <fi = ir. The dashed line close to the bot- 
tom of the figure (Pipa ~1.275xl0 -6 ) corresponds to 
the case in which only the attosecond pulses, but not the 
infrared, are acting on the system. It illustrates impres- 
sively the enormous capability of the infrared pulse to 
steer the ionized electron, thereby dramatically changing 
IPA: for 4> = the photo-association probability is ^50 
higher than without the infrared, while for <j) = 0.73-7T it 
is around 37 times smaller. The contrast between maxi- 
mum (4> = 0) and minimum (<f> — 7r) is therefore ~1850 
rendering this mechanism very efficient. Fig. [5] illustrates 
the effect of the IR field on attosecond induced photo- 
association. It shows the change in the electron trajec- 
tory due to the infrared field, for <j> = and ir, together 
with the expected trajectory of the electron under the 
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FIG. 7: Electron trajectories calculated with Eq. ((6} in the 
presence of the wells (dashed lines), the infrared field (dotted 
lines), or both of them (full lines), for carrier-envelope phases 
of <f> = or (f> = n. The position of the wells is indicated by 
the dashed lines. The origin in time is the center of the first 
attosecond pulse. 



influence of the two nuclei and the IR field. These tra- 
jectories have been calculated with Eq. where now 
Vm,(x(t)) = V + xEut(t), i.e., the interaction is given by 
the original potential V from Eq. ((2|) with the IR-ficld 
added. 

With = 0, the electron is maximally decelerated. It 
arrives later at the second well, enhancing the recombi- 
nation probability. For <f> = ir the effect is a slight ac- 
celeration: the electron arrives even earlier at the second 
well, hence the IPA probability is vanishing small. This 
can be seen even more clearly in Fig. [7] in terms of clas- 
sical trajectories of the electron for the cases shown in 
Fig.[6l From there it is clear, that the optimum situation 
for IPA (coincidence in space and time of electron trajec- 
tory, second well position and peak of second attosecond 
pulse) cannot be reached for the present parameters cho- 
sen, which was done on purpose to provide a "typical" 
situation. 



V. TRAIN OF ATTOSECOND PULSES 

The combination of two attosecond pulses and an IR 
pulse as discussed in the previous section can enhance 
the IPA probability. Yet, it is still small, especially con- 
sidering that in one dimension the recombination proba- 
bilities are unrcalistically increased due to the absence of 
wavepacket spreading in the transversal direction. One 
may overcome this problem by using more than two at- 
tosecond pulses. With more pulses in the train and a 
longer IR pulse, the transfer process is repeated each cy- 
cle, thus increasing the net probability of recombination 
at W- . From the experimental point of view, longer APT 
are even easier to generate and to handle. 

After the second pulse in the train, the population at 
the second well is non-zero. Hence, for any further pulse, 
there will be ionization from both wells. Yet, given a 
fixed ionization probability p, there is a net increase of 
population at W_, because initially the probability at 
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FIG. 8: Probability of induced photo-association at W- as a 
function of the carrier-envelope phase of the infrared pulse, 
for 2, 4 and 8 attosecond pulses in the train. Intensities: 
Iir = 10 13 W/cm 2 , and I uv = 10 15 W/cm 2 . 

W + is much larger than at . 

For n pulses, a simple calculation of the population 
probabilities at W+ and W- yields 

At = A+_ 1 {l-p) + A-^ 2 p 2 

A- = A-_ 1 (l-p)+A+_ 2 p 2 . (18) 

If initially the first well is occupied (A$ = 1) and the sec- 
ond well is empty (Aq = 0), the population of the second 
well as a function of the number of pulses n will evolve 
to a maximum whose position depends on the value of p. 
In our case, the first attosecond pulse induces an ioniza- 
tion of about 3.64% of the initial state, thus p = 0.0364 
and the maximum appears at n = 29. In Fig. [5] we show 
the IPA in W- versus the carrier-envelope phase for nu- 
merical calculation using an APT containing 2, 4 or 8 
pulses. By increasing the number of attosecond pulses, 
the absorption can be better controlled: for 2 pulses the 
enhancement factor at = with respect to the case 
without IR is 50 (as seen in Fig. [2]) increasing to 480 and 
1530 for 4 and 8 pulses, respectively. 

VI. CONCLUSIONS AND OUTLOOK 

We have introduced the process of attosecond pulse 
induced photo-association. For that purpose, a pump- 
probe scheme of attosecond pulses has been used in a 
system consisting on a dissociating H 2 molecule with the 
electron initially located in one of the ions. The pump 
pulse ionizes the electron and the probe pulse can induce 
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the reabsorption at the other ion when the delay between 
the pulses corresponds to the time the electron needs to 
reach the second ion. 

We have proven that the IPA probability can be con- 
trolled by changing the delay between the pulses. We 
have also discussed and analyzed two experimental fea- 
sible ways for enhancing this probability: the use of an 
infrared pulse to accelerate or decelerate the electron, 
and the use of a train of attosecond pulses instead of 
a single pair, to maximize the charge transfer between 
them. Both ways are shown to enhance IPA considerably. 
Furthermore, the experimental detection of reabsorption 
is facilitated through a study of the forward/backward 
asymmetry of the electronic flux along the polarization 
direction. 

Deliberately, we have chosen the laser parameters and 
the distances of, as well as the binding energies at the 
atomic centers such that the resulting dynamics can be 
described quasiclassically, providing an intuitive picture 
in terms of classical trajectories guiding the electronic 
wavepackcts. This has been important to allow for a 
straightforward interpretation of the IPA results, as pre- 
sented here. In fact, IPA dynamics is far more intricate, 
especially, outside the quasi-classical regime, e.g., if the 
excess energy of the electronic wavepacket is small (due 
to a small frequency of the attosecond pulses or a large 
binding energy, as studied for an atom in [T^L Future 
work on more than one atomic center will engage in this 
direction with the chance to uncover so far unknown ef- 
fects. 

Moreover, our ideas can be extended to more com- 
plex molecules, where several centers may exist, open- 
ing new possibilities to carry out microscopic time-of- 
flight measurements to explore microscopic spatial land- 
scapes. This will be especially useful in situations where 
the static and well developed diffraction methods fail due 
to a fast change of the microscopic landscape, e.g., a fast 
explosion or reorganization of the ions in a large molecule. 
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